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A  Nighttime  Structure  Model  of  Atmospheric  Optical 
Turbulence,  Cj ,  Derived  from  Thermosonde  and 

High  Resolution  Rawinsonde  Measurements 


1.  INTRODUCTION 


It  is  well  known  that  small  fluctuations  in  atmospheric  temperature  or  density  cause  fluctuations  in  the 
atmospheric  index  of  refraction.  These  turbulence  induced  irregularities  cause  wavefront  distortions  of  optical, 
infrared,  radio,  and  microwave  radiation.  Such  optical  distortions  are  responsible  for  random  variations  of  irradiance 
(scintillation),  image  blurring,  and  image  distortion.  Though  these  variations  are  random,  statistical  averages  can  be 
used  to  describe  their  effects.  Spatial  statistics  of  the  random  index  fluctuations  are  characterized  by  structure 
functions  and  within  the  inertial  subrange  of  isotropic  turbulence,  the  structure  function  is  proportional  to  a  power  of 
the  spatial  separation  of  the  random  variable.  The  constants  of  proportionality  for  temperature  and  refractive  index 
fluctuations  are  known  as  Cj  and  C\  respectively.  In  describing  turbulence  effects  on  optical  propagation,  is  the 

single  most  important  physical  parameter.  This  report  presents  an  altitude  dependent  thermosonde  derived  statistical 
model  for  estimating  (£  as  a  function  of  rawinsonde  measured  parameters. 


Standard  balloon  borne  rawinsondes  provide  altitude  measurements  of  temperature,  pressure,  humidity,  and 
winds  aloft.  High  resolution  rawinsondes  also  provide  a  means  for  estimating  atmospheric  temperature  gradients  and 
wind  speed  variances.  The  Phillips  Labor!  ory  tbermosonde  attached  to  a  high  resolution  rawinsonde  provides  an 

estimate  of  C?  as  a  function  of  altitude1  c£  is  evaluated  from  C2  by  the  expression  C%  =^79.9xl0-6-^-j  C2. 


Measurements  from  over  200  tbermosonde  flights  provide  a  basis  for  a  simple  quasi-empirical  rawinsonde  model.  The 
model  is  not  fundamental  in  that  neither  the  first  principles  convective  beat  fluxes  nor  eddy  diffusion  terms  are 
invoked.  Nevertheless,  the  simplicity  of  the  approach  and  the  statistical  agreement  with  direct  measurement  warrant 
consideration. 


(Received  for  Publication  25  January  1993) 


1  Brown,  J.H.,  Good,  R.E.,  Bench,  P.M.,  and  Faucher,  G.,  Sonde  Experiments  for  Comparative  Measurements  of  Optical  Turbulence, 
AFGL-TR-82-0079,  February  1982.  ADA118740 
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2.  PREVIOUS  MODELS 


2.1.  Hufnagel  Model 


Several  models  of  Cl  have  been  introduced  into  the  literature.  A  survey  by  Good  et  ai.  reviewed  the  bulk  of  these 
up  to  19882.  Two  simple  models  given  by  Hufnagel3  are: 

(2(h)- hm 

and, 

(2(h)  =  8.2  x10-5WV(A/,0°°)  +  2.7  x  10'16e^*/15°°^ 

where,  h  is  meters  above  sea  level,  and  U  is  a  root-mean-square  wind  speed  averaged  over  the  5  to  20  km  altitude 
interval.  The  second  of  these  models  is  attractive  because  of  the  RMS  wind  speed  dependency,  which  Hufnagel  states 
has  a  "good  correlation"  with  atmospheric  turbulence.  Because  these  models  are  simple,  they  do  not  provide  for  the 
thin  layering  structure  of  turbulence.  Structure  could  be  arranged  by  introduction  of  correlated  gaussian  deviates  but 
the  result  would  be  independent  of  any  particular  turbulence  profile.  Also,  the  model’s  success  in  predicting  the 
average  <2  measurement  as  determined  by  thermosondes  often  proves  unsatisfactory2. 

22.  NOAA  Model 


A  probabilistic  model  developed  by  the  Aeronomy  Laboratory  at  NOAA  gives  <2  in  terms  of  a  multivariate 
probability  density  function  (pdf)4-5'6.  The  independent  parameters  are  N,  S,  q ,  and  L,  which  are  respectively  the 
Brunt-Vaisala  frequency,  vertical  shear  of  the  horizontal  winds,  vertical  gradient  of  the  specific  humidity,  and  the 
outer  scale  of  the  inertial  subrange.  The  reduced  model  for  a  dry  atmosphere  is  given  by: 

C2(dry)  =  dSpsf_^C  dN2pN(NZf 


2Good,  R.R,  et  al.  (1988)  Atmospheric  Modeb  of  Optical  Turbulence,  SPIE  V928  Modeling  of  the  Atmosphere,  Rothman,  L.S.  ed.,  pp  165- 

186. 


3Wolfe,  W.L.  and  Zississ,  G.J.,  editors,  ’The  Infrared  Handbook,  Chapter  6’,  Hufnagel,  R.E.,  amh.,  Propagation  Through  Atmospheric 
Turbulence ,  pg.  6-14,  Department  of  the  Navy,  1978. 

4Warnock,  J.M.  and  VanZandt,  T.E.,  (1985)  A  Statistical  Model  to  Estimate  The  Refractmty  Turbulence  Structure  Constant  C%  In  Free 
Atmosphere,  NOAA  technical  memorandum  ERL  AL-10,  Aeronomy  Laboratory,  NOAA,  Environmental  Research  Laboratories. 

-’VanZandt,  T.  E.,  Gage,  K.  S.,  and  Warnock,  J.  M.,  (1981).,  An  Improved  Model  for  The  Calculation  of  profiles  of  C%  and  e  In  The  Free 
Atmosphere  From  Background  Profiles  Of  Wind,  Temperature,  And  Humidity,  Proc.  20th  Conf.  on  Radar  Meteorology,  Nov.  30-Dec  3, 
Boston,  MA,  129-261. 

^VanZandt,  etaL,  Vertical  Profiles  Of  Refractivity  Turbulence  Structure  Constant:  Comparison  Of  Observations  By  The  Sunset  Radar  With 
A  New  Model,  Radio  Science,  13,:819-829. 
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pressure,  T 


temperature,  and  where 


where,  C,=2.8,  Af0  =  C,— ,  02  2-77.6x10^,  P  = 

gT 

pL,  ps,  and  pN ,  are  respectively  the  probability  distributions  of  L,  S,  and  N. 

At  a  smoothed  resolution  of  500m,  the  NOAA  model  has  been  compared7  very  favorably  to  radar  and 
thermosonde  measurements  of  but  its  application  remains  quite  complicated.  Also  the  standard  deviations  of  the 
distributions  are  determined  from  certain  empirical  scaling  factors,  some  of  which  depend  upon  location  and  perhaps 
measurement  resolution.  In  particular,  the  constant  in  the  equation  giving  the  distribution  of  wind  shears  depends  upon 
"the  range  of  scales  important  in  the  onset  of  turbulence."8  While  recommended  for  detailed  radar  and  other  low 
resolution  measurements,  this  technique  requires  further  evidence  for  its  applicability  toward  producing  high 
resolution  C2  structure. 

23.  Dewan  et  al.  Model 


Another  statistical  model  developed  by  Dewan,  et  al.2,9  expresses  in  terms  of  the  well  know  Tatarski 

relationship10  cf  =  2.8Af2z/^,  where  the  factor  M  gives  the  gradient  of  the  index  of  refraction,  M  =  -79  x  KF6  . 

gT 

Dewan  et  al.  estimate  L  through  a  regression  model  of  the  form:  L£>g1o[l(z)]  =  -l  +  — Y(z)  where,  y(z)  =  Cj  +  (^5, 

and  where  S  is  the  shear  evaluated  at  the  standard  rawinsonde  altitude  resolution  for  wind  speed  of  300m.  Very  high 
resolution  smoke  frail  wind  shear  data  (10m)  was  employed  in  the  regression  analysis  to  evaluate  the  constants  Cj  and 
C2.  In  this  case,  a  layer  thickness  was  assigned  to  each  region  where  the  smoke  frail  microshears  exceeded  a  critical 
value  (that  is  the  shear  necessary  to  induce  turbulence  based  on  the  Richardson  number  equal  to  1/4  criterion).  In  turn 
these  layer  thicknesses  were  related  to  300m  resolution  rawinsonde  shears.  Results  from  this  model  will  be  compared 
to  a  thermosonde  deduced  model  which  is  developed  in  the  following  section.  A  reduced  form  of  the  Dewan  et  al. 
model  takes  the  form:  L(z)  =  1.5e5u5S^  in  the  troposphere,  and  L(z)  =  0.24e63,925^  in  the  stratosphere. 

2  A.  Thermosonde  Model 


The  thermosonde  model  is  developed  from  the  Tatarski  expression,  c£  =  2.8M2L%,  in  the  manner  of  the 
previous  model,  except  that  the  variables  governing  M  are  measured  in-situ  as  a  function  of  altitude.  That  is, 

C2(z)  =  2.8[Af(z)]2[L(z)j3.  In  this  expression,  M(z)  =  -79xl 0-6  where  the  altitude  dependent 

parameters  are:  P(z)  (ambient  pressure  in  mb),  T(z)  (ambient  temperature  in  Kelvin),  and  coB(z)  (local  Brunt-Vaisala 


7  7 

'War nock.  J.M.,  etal.,  Comparison  Among  Clear-Air  Radar,  Thermosonde  and  Optical  Measurements  and  Model  Estimates  of C„  Made  in 

Very  Flat  Terrain  Over  Illinois,  Middle  Atmospheric  Program  Handbook,  Liu,  C.H.,  ed.,  V28,  June.1989 

8Green,  J.L.,  et.  al.  Comparisons  of  Refractivity  Turbulence  Estimates  from  the  Flatland  VHF  Radar  with  Other  Measurement  Techniques, 
AMS,  24th  Conf.  on  Radar  Meterol.,  Tallahasee  FL,  March  1989. 

9  Dew  an,  E.M.,  Good,  R.E.,  Beland,  R.,  and  Brown,  J.H.,  "A  Model  for  C2  (Optical  Turbulence)  Profiles  from  Radiosonde  Data".  PL- 
TR-93-2043 

10Tatarski,  V.  I.,  (1961)  Wave  Propagation  in  a  Turbulent  Medium,  McGraw-Hill 


3 


frequency).  The  square  of  the  Brunt- Vaisala  frequency  is  estimated  from  in-situ  temperature  measurements  by  the 

expression  tol(z)  =  — 7-rf  — ' —  ■ +  r )  where  r  is  the  dry  adiabatic  lapse  rate,  0.0098  K/m.  Given  these  quantities,  L 
T{z){  dz  ) 

was  estimated  from  the  thermosonde  measurements  by  the  expression: 


Lfhermosonde  ~ 


y* 


Cf  (smooth) 


2.8  x 


8 


.  A  sum  of  the  squares  minimization  between  i-thennosonde  and  over  tw0 


hundred  different  combinations  of  rawinsonde  local  variables  and  derived  variables  was  employed  to  find  a  regression 
form  that  yielded  the  smallest  residual  error.  The  variables  of  most  significance  resulted  from  a  combination  of  the 
local  values  of  the  Brunt-Vaisala  frequency,  wind  speed  variance,  and  their  averages.  In  particular,  the  measurements 
led  to  two  suitable  empirical  models  given  by  the  following: 


Model  (1) 

^og io(*Mdi)  =  °i +  <hLogl0(°l) + OjLo&oCMb  +  aA ) 
Model  (2) 


**&oO«.i2>=*k  +b2te8l0(£)+b3Logl0(!bd^) 

av  0>8 

where  a J  is  an  estimate  of  the  local  wind  speed  variance  and  where  the  overbars  denote  their  averages  over  the  entire 

troposphere  or  stratosphere,  depending  upon  the  data  region  {i.e.  different  sets  of  coefficients  were  determined  for  the 
two  different  regions).  The  estimated  quantities  of  to|  and  aj  were  derived  from  high  resolution  rawinsondes  by 

treating  the  data  over  150m  altitude  intervals.  Linear  interpolation  of  the  20m  raw  data  was  employed  whenever  it  was 
necessary  to  "line  up"  values  at  the  same  altitude  levels.  Special  consideration  was  given  toward  estimating  cog  and 
Gy .  In  calculating  Gig,  the  smoothed  estimate  of  the  temperature  gradient,  dT/dz,  was  found  by  performing  a  five 
point  running  average  over  the  20  m  high  resolution  rawinsonde  temperature  and  then  performing  a  least  squares 
parabolic  fit  over  13  consecutive  values  of  the  resulting  array.  The  temperature  derivative  was  found  from  the  slope  of 
the  estimated  parabola  near  the  seventh  value  used  in  finding  the  fit.  In  calculating  Gy,  a  running  average  of  the 
LORAN-derived  wind  speed  over  150  m  was  subtracted  from  the  local  values.  Then  Gy  was  estimated  by  triangular 
weighting  the  sum  of  the  squares  of  the  result,  that  is. 
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X(v-W2  *  Weighty; 

Oviert.  =  ove,'150m.— - - - .  Since  the  Govind1 1  filtering  process  employed  to  calculate  wind  speed 

Weightfnang 
over  150m 

from  the  transmitted  LORAN  time  delays  produced  altitude  resolutions  much  thicker  than  20m,  some  of  the  low  signal 
contribution  to  a2  is  noise  induced.  There  appears  to  exist,  however,  sufficient  strong  velocity  fluctuation 
contrib”tions  to  c2  and  evident  tracking  between  major  peaks  in  a2  and  to  warrant  use  of  this  parameterization. 
Furthermore,  use  of  a2  in  the  thermosonde  model  reduces  the  residual  error  between  measurement  and  model. 


3.  ANALYSIS 


Fifteen  nighttime  thermosonde  flights  were  used  to  calculate  the  model  constants.  These  were  launched  from 
Pennsylvania  State  University  in  May  1986.  A  typical  flight  (Launch  L4007)  is  characterized  in  Figure  1.  The  graph 
shows  three  panels,  where  the  left  panel  contains  plots  of  ambient  temperature  and  humidity,  the  right  panel  contains 
plots  of  wind  speed  and  direction,  and  the  center  panel  contains  a  plot  of  high  resolution  C^.  Note  that  the  horizontal 


dashed  lines  indicate  the  altitudes  of  boundary  layer  temperature  and  humidity  inversions  and,  higher  up,  the 
tropopause.  These  levels  define  the  troposphere  and  stratosphere  for  the  purposes  of  this  report  Overdrawn  on  the 


plot  is  a  solid  smooth  curve  that  shows  a  simple  empirical  model  for  the  entire  Pennsylvania  State  University 

\4t\' 

nighttime  data  set.  This  model  takes  the  form:  Logi0[c£ (z)]  -a+bz+cz2  +dv- 


It  will  be  called  the  Penn 


State  Empirical  Model  (PSEM)  and  it  will  be  used  as  a  visual  reference  in  comparisons  with  Model  (1)  and  Model  (2). 
The  next  series  of  figures  (Figures  2-5)  depict  a  case  in  the  development  of  Model  (1).  The  left  panel  in  Figure  2 
shows  the  L4007  raw  (20m)  altitude  profile  of  C2  while  the  middle  and  right  panels  compare  the  smoothed  (150m) 
C2  with  C2  derived  from  Model  (1).  As  should  be  expected,  the  model  is  smoother  than  the  thermosonde  data  but 


quantitatively  and  structurally  the  agreement  is  very  good.  Figure  3  compares  L  derived  from  the  thermosonde  data 
(left  panel)  with  Model  (1)  (right  panel).  Alignment  of  the  peaks  is  reasonable.  Figure  4  shows  the  altitude  profiles  of 
the  measured  rawinsonde  parameters  (0B  (left  panel)  and  cv  (right  panel).  Fifteen  similar  profiles  of  (Dg,  a2,  and  cf 
were  used  to  evaluate  the  coefficients  of  Model  (1)  and  Model  (2).  The  middle  panel  of  Figure  4  shows  the  model 
profile  for  L4007.  This  may  be  compared  to  the  PSEM  depicted  as  the  smooth  curve  overlay  and  to  the  raw  of 


Figure  1 .  Again  there  is  good  quantitative  agreement  between  model  and  the  smoothed  measurement  A  comparison  of 
L  derived  for  Model  (1)  and  the  L4007  thermosonde  derived  measurement  of  L  is  shown  in  the  scatter  plot  of  figure  5. 
The  left  and  right  panels  are  for  tropospheric  and  stratospheric  data  respectively.  A  45  degree  line  on  these  plots 
indicates  exact  agreement.  The  range  of  variation  of  L  is  shown  to  be  wider  in  the  troposphere  than  in  the  stratosphere, 
nevertheless,  the  data  tends  to  fall  about  the  45  degree  slope.  Other  individual  cases  are  presented  later  in  this  report 
for  comparison.  Figure  6  presents  a  comparison  of  L  derived  for  Model  (1)  and  all  fifteen  thermosonde  profile  derived 
measurements  of  L.  Figure  7  does  the  same  for  Model  (2).  Here  the  x's  depict  the  median  of  the  binned  data  and  the 
error  bars  depict  the  standard  deviations.  Numbers  above  the  error  bars  indicate  the  number  of  points  in  the  bin  and 


*  ^Govind,  P.K.,  (1975)  Omega  Windfinding  Systems,  J.  AppL  Meteor .,  V14: 1503-1515 


5 


numbers  above  the  plots  indicate  bow  many  points  lay  outside  three  standard  deviations.  Again  the  data  tend  to  fall 
about  the  45  degree  slope. 

As  evaluated  from  the  regression  analysis,  the  constants  and  their  standard  deviations  for  Model  (1)  and  Model 
(2)  are  given  in  Table  1 . 


Table  1.  Regression 


Constants  for  Model  (1)  and  Model  (2) 

Model  (1) 

Model  (2) 

Troposphere 

Stratosphere 

Troposphere 

Stratosphere 

al 

-4.004 

-2.68 

bl 

-0.1684 

-0.276 

aal 

0.0457 

0.0592 

°bl 

0.00489 

0.00582 

a2 

0.1729 

0.2285 

*2 

0.1550 

0.2136 

®a2 

0.0146 

0.0168 

°b2 

0.0148 

0.0173 

a3 

-2.00 

-1.473 

^3 

-2.072 

-1.385 

°a3 

0.0225 

0.0356 

°b3 

0.0236 

0.0337 

34 

-0.00054 

-0.00295 

b4 

-0.00054 

-.00307 

CTa4 

.0000041 

.000026 

°b4 

.0000036 

.0000181 

_  (o}P 

L  may  be  written,  for  example,  for  Model  (1),  in  the  troposphere  as:  L  =  10  - - which  suggests 

(a>*-.00054)z 

that  L  may  scale  with  the  variance  of  the  wind  speed  and  inversely  with  the  Brunt-Vaisala  frequency.  In  fact  the 

I  _i 

outerscale  of  inertial  range  turbulence  Lq  is  propontional  to  e2cob2,  where  e  is  the  turbulence  dissipation  rate  and 
where  e2  is  proportional  to  ,  where  v'2  is  the  variance  of  the  small  scale  wind  speed  fluctuations12’13,  so  that 


12Weiiistock,  (1978)  Vertical  Turbulent  Diffusion  in  a  Stably  Stratified  Fluid,  J.  Atm.  Sd.,  V35: 1022. 

'^Hocking,  W.K.,(1985)  Measurement  of  turbulent  energy  dissipation  rates  in  the  middle  atmosphere  by  radar  techniques:  A  review. 
Radio  Science.  20.,  (No.  6),  1403-1422. 
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-T 

Zo  - - .  Since  the  rawinsonde  measures  large  scale  wind  speed  and  not  v,z ,  Model  (1)  and  Model  (2)  estimate  the 

CO® 

microscale  through  the  scaling  exponents  of  and  to*. 


4.  OTHER  CASES 

4.1.  Qualitative  Comparison 


Figures  1-5  presented  in  the  preceding  section  was  a  Pennsylvania  State  case  for  Model  (1).  This  section  presents 
additional  cases  for  Pennsylvania  State  and  for  other  locations  and  includes  cases  for  both  Model  (1)  and  Model  (2). 
Figures  8-11  show  the  same  case  as  before  (L4007)  but  for  Model  (2).  There  exist  only  minor  differences  between  the 
c£  results  for  Model  (1)  and  Model  (2)  but  differences  in  the  magnitude  of  L  are  evident  The  intended  purpose  of 
introducing  the  average  parameters  of  a>|  and  ctj  in  Model  (2)  is  to  provide  an  adjustment  for  low  winds  aloft  versus 
jet  stream  conditions.  Figures  12-19  show  another  Pennsylvania  State  University  thermosonde  flight  (L4044).  Model 
(1)  is  presented  in  Figures  13-16  while  Model  (2)  is  presented  in  Figures  17-19.  This  flight  differs  from  L4007  in  the 
very  low  C%  values  in  the  troposphere  from  6-10  km  and  also  in  the  peak  enhancement  of  Cz  just  above  the 
tropopause  at  13  km.  Figure  15  shows  that  the  model  tracks  the  low  tropospheric  values  of  as  well  as  the  13  km 

enhancement.  Examination  of  the  Brunt-Vaisala  frequency  plot  and  the  wind  speed  variance  plot  in  Figure  15  reveals 
that  the  enhancement  can  be  traced  to  the  enhancement  of  these  parameters  at  13  km. 

Figures  20-28  show  a  thermosonde  flight  (L1014)  conducted  at  Champaign,  Illinois  on  the  night  of  14  June  , 
1988.  Figures  21-24  show  the  results  for  Model  (1).  Figure  25  presents  a  plot  of  the  shear  and  Richardson  number  for 
comparison  to  the  model.  Figures  26-28  show  the  results  for  Model  (2).  Even  though  the  model  constants  are  those 
that  were  evaluated  for  the  Pennsylvania  State  campaign,  visual  examination  will  reveal  very  good  agreement  between 
the  thermosonde  and  the  model  C^.  Note  particularly  the  level  of  agreement  between  L  of  the  data  and  L  of  the 

model  in  Figure  22,  especially  at  7  km.  Again  this  agreement  can  be  traced  to  enhancements  in  the  Brunt-Vaisala  fre¬ 
quency  and  wind  speed  variance.  Although  the  Richardson  number  criterion  for  turbulence  is  not  invoked  in  the 
present  model,  there  does  exist  a  strong  feature  in  the  Richardson  number  plot  in  Figure  25  at  7  km.  This  observation 
suggests  that  the  relationship  of  local  microshear  and  Brunt-Vaisala  frequency  can  be  represented  by  scaled 
rawinsonde  measurements  of  Brunt-Vaisala  frequency  and  wind  speed  variance  as  depicted  in  Figure  23. 

Figures  29-33  show  a  thermosonde  flight  (M0516)  conducted  at  a  desert  site  in  New  Mexico  on  the  night  of  5 
September  ,  1984.  Figures  30-32  show  the  results  for  Model  (1)  and  Figure  33  shows  the  result  for  Model  (2). 

Again  the  same  model  constants  are  used  but  applied  to  a  climatologically  different  site.  Comparison  of  Figures  29 
with  figures  32  and  33  show  good  agreement.  Model  (2)  shows  slight  improvement  in  the  stratosphere. 

Figures  34-37  show  another  thermosonde  flight  (M6385)  conducted  at  the  same  desert  site  as  M0516  but  at  a 
different  time  of  year  (3  March  1985).  Figures  35-36  show  the  results  for  Model  (1)  and  Figure  37  shows  the  result 

for  Model  (2).  This  case  presents  a  challenge  to  the  model  because  no  winds  aloft  measurements  were  available. 
Lacking  the  measured  wind  speeds,  the  model  was  calculated  for  an  arbitrary  RMS  wind  speed  of  1  m/s.  An 
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examination  of  the  agreement  between  the  thennosonde  profile  in  Figure  34  and  the  model  in  Figure  36  suggests 
that  the  scaled  Brunt-Vaisaia  frequency  alone  may  be  sufficient  to  characterize  Q*  in  the  absence  of  very  strong  shear 
layers.  This  concept  will  be  examined  further  in  the  next  section. 

4.2.  Quantitative  Comparison 


A  parameter  that  is  sometimes  used  in  describing  optical  turbulence  is  the  isoplanatic  angle.  This  is  an  altitude 
weighted  integral  of  that  has  significance  in  optical  imaging  but  which  can  also  save  as  a  basis  of  quantitative 

comparison  between  measurements  of  different  instruments  or  between  measurement  and  model.  The  isoplanatic 
angle  is  evaluated  as : 


where  h  is  height  above  ground  and  k  is  wavenumber.  For  the  purposes  of  this  section,  the  limits  of  integration  are 
taken  over  the  measured  thennosonde  altitude  region.  Tables  2-3  list  60  for  the  thermosonde  measurements  and 
models  at  various  locations. 


Table  2.  Isoplantic  Angle  Calculations  for  the  Pennsylvania  State 
University  Campaign.  Comparison  of  Measurement  to  Model 


Pennsylvania  State  University  Campaign 

60X  10*  rad 

Flight# 

Measurement 

Model  (1) 

Model  (2) 

L4007 

6.6 

7.9 

7.7 

L4012 

4.7 

8.2 

7.4 

L4014 

6.1 

6.4 

8.5 

L4018 

8.0 

8.3 

8.2 

L4019 

8.5 

9.1 

8.8 

L4029 

5.5 

7.4 

7.3 

L4031 

6.4 

8.4 

8.5 

L4032 

4.9 

7.6 

7.7 
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Table  3  Isoplanatic  Angle  Calculations  for  Champaign  Illinois 
Campaign.  Comparison  of  Measurement  to  Model 


Champaign,  Illinois  Campaign 

0OX  10*  rad 

Flight# 

Measurement 

Model  (1) 

Model  (2) 

L0553 

10. 

8.1 

8.7 

II 

12. 

9.3 

9.3 

L1006 

11. 

9.1 

8.3 

L1014 

8.0 

8.5 

8.6 

L3990 

6.8 

9.9 

9.9 

Average 

9.6 

$.0 

9.0 

Std.  Dev. 

1.9 

.62 

.57 
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Table  4.  Isoplanatic  Angle  Calculations  for  New  Mexico  Campaign 
-  Sept  1984.  Comparison  of  Measurement  to  Model 


New  Mexico  Campaign  •  Sept.  1984 

e0X  10*  rad 

Measurement 

M0516 

11. 

8.0 

9.0 

M0522 

7.5 

6.8 

6.1 

M1880 

11. 

9.3 

10. 

M1886 

3.5 

7.9 

8.9 

Average 

8.2 

8.0 

8.6 

Std.  Dev. 

3.1 

.87 

1.5 

Table  5.  Isoplanatic  Angle  Calculations  for  New  Mexico  Campaign 
-  March  198S.  Comparison  of  Measurement  to  Model 


New  Mexico  Campaign  -  March  1985 

e„X  10*  rad 

■sph 

Measurement 

Model  (1) 

M6381 

7.2 

12. 

8.8 

M6385 

12. 

12. 

8.7 

M6386 

6.3 

15. 

9.0 

M6408 

5.1 

8.9 

8.5 

M6417 

3.9 

9.6 

8.6 

M6563 

8.3 

11. 

10. 

Average 

7.1 

11. 

9.0 

Std.  Dev. 

2.5 

1.9 

23 
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The  Pennsylvania  State  data  of  Table  2  indicate  that,  on  average,  both  models  overestimate  0j  by  about  20 
percent  and  reduce  the  variance  of  ©0  somewhat  Model  (1)  however  provides  more  variance  than  Model  (2).  The 
Champaign,  Illinois  data  of  Table  3  indicate  that,  on  average,  both  models  underestimate  60  by  about  6  percent  Model 
(1)  of  Table  4  underestimates  60  by  3  percent  while  Model  (2)  overestimates  by  4  percent.  Here,  Model  (2)  provides 
more  variance.  The  desert  site  March  1985  data  of  Table  5  display  the  most  error.  Model  (1)  overestimates  ©0  by  61 
percent  while  Model  (2)  overestimates  60  by  27  percent.  The  problem  seems  to  rest  in  the  fact  that  no  winds  aloft  data 
was  used  for  this  particular  series.  In  order  to  apply  the  model,  the  RMS  wind  speed  for  these  flights  was  arbitrarily  set 
to  1  m/s  but  the  poor  result  indicates  that  al  should  be  included  after  all. 

43.  Comparison  with  the  Dewan  et  al  Model 

As  described  earlier,  the  Dewan  et  al.  model  is  a  scaled  wind  shear  model.  In  this  section  we  will  examine  the 
Dewan  et  al.  model  in  toms  of  the  Pennsylvania  State  University  campaign  data.  Since  the  model  coefficients  were 
derived  for  300  m  resolution,  the  applicable  rawinsonde  parameters  of  temperature,  Brunt-Vaisala  frequency,  and 
wind  shear  were  smoothed  to  300  m  resolution.  Figure  38  shows  the  wind  shear  profile  of  flight  L4007  that,  when 
applied,  results  in  the  model  plots  of  Figures  39-42.  Examination  of  Figures  39  and  40  reveals  very  little  tracking  of 
the  small  scale  profile  structure  between  the  thermosonde  c}  and  model  C f  or  between  the  thermosonde  smoothed 

estimate  of  L  and  L  estimated  from  the  model.  Figure  42  also  indicates  poor  tracking  of  the  thermosonde  estimate  of  L 
and  L  estimated  from  the  model.  On  the  other  hand,  examination  of  Figure  41  reveals  that  the  average  profile  agrees 
well  with  the  flight  data.  Also  some  of  the  large  scale  structure  (2-5  km)  tracks  the  thermosonde  profile  fairly  well. 

Another  flight  (L4012)  is  examined  in  Figures  43-47.  Lack  of  tracking  is  again  evident  in  Figures  44, 45,  and  47. 
It  would  appear  that  the  Brunt-Vaisala  term  accounts  for  the  large  scale  tracking  of  in  Figures  41  and  46. 

0O  was  also  calculated  for  the  Dewan  et  al.  me  lei  and  the  results  are  tabulated  in  Table  6.  On  average,  the  model 
underestimates  0O  by  32  percent.  Considering  the  single  parameter  estimate  of  the  microscale  shear,  the  model 
performs  well  in  estimating  the  C%  average  altitude  profile,  fairly  in  estimating  0O,  but  poorly  in  tracking  small  scale 
structure  in  the  Q  profile.  It  appears,  therefore,  that  the  microscale  shear  estimate  is  insufficient  for  predicting  Q 
structure  with  thicknesses  less  than  2  km. 
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Table  6.  Isoplan&tic  Angie  Calculations  far 
Dew  an  et  al.  Model  Comparison  of  Measurement 
to  Model 


80  X  10"6  rad 

Flight# 

Measurement 

Dewan  et  al. 
Model 

L4007 

6.6 

4.3 

L4012 

4.7 

2.7 

L4014 

6.1 

.25 

L4018 

8.0 

5.6 

L4019 

8.5 

mem 

L4029 

5.5 

3.2 

L4031 

6.4 

4.8 

L4032 

4.9 

4.8 

L4033 

6.0 

5.4 

L4035 

4.7 

3.9 

L4037 

9.0 

5.4 

L4042 

9.0 

4.7 

L4043 

5.1 

4.1 

L4044 

6.9 

5.6 

L4045 

7.9 

3.2 

without  L4014 

Average 

6.6 

4.5 

Std.  Dev. 

1.5 

.97 
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5.  CONCLUSIONS 


A  comprehensive  study  of  thennosonde  data  resulted  in  a  simple  nighttime  quasi-empirical  rawinsonde  model 
that  takes  the  form: 

Q2(z)  =  2.8(Af(z)f[^z)F, 
where 


M(z)  =  -79  xlO~* 

gT(z) 


and  where  L(z)  =  C\ 


The  model  for  L(z)  is  an  exponentially  scaled  estimate  of  the  theoretical  description  and  it  uses  high  resolution  rawin¬ 
sonde  data  to  estimate  the  altitude  dependent  parameters  c2  and  to*.  The  c£  model  is  independent  of  location  but 

applicable  to  a  variety  a  different  climatological  conditions.  Compared  to  other  models,  the  thennosonde  model  is 
either  much  less  complex  or  shows  better  fidelity  to  thermosonde  derived  measurements. 

Thermosondes  are  reliable  instruments  for  obtaining  tropospheric  and  stratospheric  optical  turbulence  data. 
When  optical  turbulence  is  important  in  characterizing  new  sites  for  vertical  paths,  the  described  model  may  be 
employed  to  estimate  c£.  Many  low  cost,  high  resolution  rawinsondes  can  be  deployed  rapidly  to  construct  a 
turbulence  climatology.  Where  high  resolution  data  already  exists,  the  proposed  model  offers  significant  time, 
robustness,  and  cost  advantages  in  exchange  for  a  small  computational  burden. 
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Figure  2.  Thetmosonde  cf  Measurement  (Raw  and  Smoothed)  Profiles  Compared  to  Model  (1)  C*  Profile.  Flight  L4007 
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Figure  4.  Brunt- Vaisala  Frequency  and  RMS  Wind  Speed  Profiles  Derived  From  Measurement  Compared  With  C*  Profile  from  Model  (1).  Flight  L4007 
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Figure  6.  Binned  Scatter  Plot  of  Data  from  Entire  Pennsylvania  State  University  Campaign.  "L"  for  smoothed  thermosonde  measurements  compared  with 
“L"  from  Model  (1).  A  45  degree  slope  represents  perfect  agreement.  Error  bars  represent  the  standard  deviation  of  the  data  in  each  bin.  Numbers 
above  the  plots  are  the  number  of  points  falling  outside  three  standard  deviations.  Left-hand  plot  is  troposphere.  Right-hand  plot  is  stratosphere 
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Figure  7.  Binned  Scatter  Plot  of  Data  From  Entire  Pennsylvania  State  University  Campaign.  "L"  for  smoothed  thermosonde  measurements  compared  with 
"L"  from  Model  (2).  A  45  degree  slope  represents  perfect  agreement.  Error  bars  represent  the  standard  deviation  of  the  data  in  each  bin.  Numbers 
above  the  plots  are  the  number  of  points  falling  outside  three  standard  deviations.  Left-hand  plot  is  troposphere.  Right-hand  plot  is  stratosphere 
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Figure  8.  Model  (2)  C,2  Profile.  Flight  L4007 
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Figure  1 1 .  Scatter  Plot  of  "L*  For  Smoothed  Thermosonde  Measurements  Compared  with  "L"  From  Model  (2).  Flight  14007.  Leftmost  plot  for 
Rightmost  plot  for  stratosphere 


but  for  Pennsylvania  State  University  Flight  L4044 
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1 3.  Same  as  Figure  2  but  for  Pennsylvania  State  University  Flight  L4044 


Figure  14.  Same  as  Figure  3  but  for  Pennsylvania  State  University  Flight  L4044 


figure  15.  Same  as  Figure  4  but  for  Pennsylvania  State  University  Flight  L4044 


1 7.  Same  as  Figure  9  but  for  Pennsylvania  State  University  Flight  L4044 
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Figure  19.  Same  as  Figure  1 1  but  for  Pennsylvania  State  University  Flight  L4044 
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Figure  22.  Same  as  Figure  3  but  for  Champaign,  Illinois  Flight  L1014,  June  1988 


23.  Same  as  Figure  4  but  for  Champaign,  Illinois  Flight  L1014,  June  1988 


Figure  24.  Same  as  Figure  S  but  for  Champaign,  Illinois  Flight  L1014,  June  1988 
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Shear  and  Richardson  Number  Profiles  Calculated  for  Champaign,  Illinois  Plight  L1014,  June  1988 
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Figure  26.  Same  as  Figure  8  but  for  Champaign,  Illinois  Flight  LI0I4,  June  1988 
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Figure  27.  Same 
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Figure  29.  Same  as  Figure  1  but  for  New  Mexico  Flight  M0516,  September  1984 
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Figure  30.  Same  as  Figure  2  but  for  New  Mexico  Flight  M05 1 6,  September  1 984 


but  for  New  Mexico  Flight  M0516,  September  1984 
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Figure  32.  Same  as  Figure  4  but  for  New  Mexico  Flight  M05 16,  September  1984 


10  but  for  New  Mexico  Flight  M0516,  September  1984 
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Figure  34.  Same  as  Figure  1  but  for  New  Mexico  Flight  M6385,  March  1985 


Figure  35.  Same  as  Figure  2  but  for  New  Mexico  Flight  M6385,  March  1985 
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Figure  36.  Same  as  Figure  4  but  for  New  Mexico  Flight  M6385,  March  1985. 


o  in  o  in  o  in  o 

m  CM  n  —  — 


nsw  wm)  aaruum 


50 


Figure  37.  Sarny  as  Figure  10  but  for  New  Mexico  Flight  M6385,  March  1985 
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Figure  38.  Smoothed  300  m  Rawinsonde  Shear  for  Pennsylvania  State  University  Flight  L4044 
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Figure  46.  Same  as  Figure  4 1  but  for  Pennsylvania  State  University  Flight  L4012,  May  1986 
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le  as  Figure  42  but  for  Pennsylvania  State  University  Flight  L4012,  May  1986 
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